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ABSTRACT 

We study the hnear evolution of small perturbations in self-gravitating fluid 
systems with magnetic fields. We consider wave-like perturbations to nonuniform 
filamentary and sheet-like hydrostatic equilibria in the presence of a uniform par- 
allel magnetic field. Motivated by observations of molecular clouds that suggest 
substantial nonthermal (turbulent) pressure, we adopt equations of state that are 
softer than isothermal. We numerically determine the dispersion relation and the 
form of the perturbations in the regime of instability. The form of the dispersion 
relation is the same for all equations of state considered, for all magnetic field 
strengths, and for both geometries examined. We demonstrate the existence of 
a fastest growing mode for the system and study how its characteristics depend 
on the amount of turbulence and the strength of the magnetic field. Gener- 
ally, turbulence tends to increase the rate and the length scale of fragmentation. 
While tending to slow the fragmentation, the magnetic field has little effect on 
the fragmentation length scale until reaching some threshold, above which the 
length scale decreases significantly. Finally, we discuss the implications of these 
results for star formation in molecular clouds. 
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1. Introduction 



It is widely believed that most current galactic star formation occurs in molecular clouds 
(cf. [Shu, Adams, fc Lizano 1987] for a review). Observations of these clouds reveal rich and 
complex structure (e.g., Myers 1991 ; Blitz 1993 ). In particular, filamentary and sheet-like 
structures are common. These structures exhibit a large degree of substructure of their 
own (e.g., Schneider fc Elmegreen 1979| ; |de Geus, Bronfman, fc Thaddeus 19901 ; |Houlahan| 
Scalo 1992 ; Wiseman fc Adams 1994 ; Wiseman 1996 ). Clumps and cores that provide 
locations for new stars to form are often found along filaments and sheets. Thus, the study 
of the evolution of molecular clouds and their substructure is important for understanding 
the processes involved in present day star formation. The goal of this paper is to further our 
understanding of molecular cloud substructure by studying gravitational instabilities in the 
presence of both turbulence and magnetic fields. 

Most previous work has focused on the case of wave motions and instabilities in uniform 
density fluids (from [Jeans 19"2^ to Pewer 1970| , [Langer 1978| , Pudritz 1990|) . Some previous 
work on wave motions in isothermal molecular cloud filaments and sheets has been done; this 
work begins with clouds in hydrostatic equilibrium and hence nonuniform density. Many ob- 
servations show clumps which appear nearly equally spaced along filaments (e.g., [Schneider 
fc Elmegreen 1975| ; [McBreen et al. 1979| ; putrey et al. 199TD ; this finding has led to the idea 
that the clumps may arise from a gravitational instability with a particular length scale. 
Less frequently, it has been proposed that the clumps may be peaks of density waves which 
propagate along the filament. In the linear regime, both cases are treated by a linear pertur- 
bation analysis. Larson (1985) gives a good review of the early progress in doing this type of 
analysis. More recent work includes Nagasawa (1987), who performed such calculations for 
an idealized isothermal filament with an axial magnetic field (see also jMiyama, Narita,~^ 



Hayashi 1987 for a discussion of sheets). Others ( Nakamura, Hanawa, fc JNakano 1991 , 1993; 
[Matsumoto, JNakamura, fc Hanawa 1994] ) have added further embellishments to this model 
(e.g., rotation). Until recently, however, these studies have typically neglected the effects of 
turbulence and non-isothermal equations of state. Gehman et al. (1996; hereafter Paper 1[ ) 
introduced a turbulent equation of state to the problem, but did not include the effects of a 
magnetic field. When including a magnetic field, one is generally free to choose the form of 
the unperturbed magnetic field. The results may depend strongly on the form assumed. In 
this present work, we consider a general barotropic equation of state of the form p = P{p) 
along with a uniform magnetic field. We concentrate this present discussion on equations of 
state which are softer than isothermal, since these may be the most relevant for molecular 
clouds (see ^0). 



A secondary motivation for this work is to generalize the current theory of star forma- 
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tion, which typically begins with spherically symmetric clouds (e.g., Larson 1972 ; Shu 1977 
lerebey, Shu, fc Cassen 1984| ). Models of star formation built upon these collapse calcula- 



tions are reasonably successful and predict spectral energy distributions of forming stars that 
are in agreement with observed protostellar candidates ( Adams, Lada, fc Shu 198"^ ; Butner et 



al. 1991| ; [Kenyon, Calvet, fc Hartmann 199"3| ). However, departures from spherical symmetry 
are expected to occur on larger spatial scales, and the collapse of isothermal self-gravitating 
sheets ( [Hartmann et al. 1994|) and filaments ([Inutsuka fc Miyama 1993| ; [Nakamura, Hanawa 



Nakano 1995|) has begun to be studied. The calculations of this paper help to determine 
the length scales of fragmentation of molecular cloud sheets and filaments, and the initial 
conditions for protostellar collapse in nonspherical molecular cloud structures. 

This paper is organized as follows. We present the general mathematical framework 
and discuss equations of state in The equilibrium solutions of interest are presented and 
discussed in §|[ In §^ we perform a general linear perturbation analysis, the solutions of 
which are presented and examined in §^ We conclude in with a discussion and summary 
of our results. 



2. General Formulation 

We present here the basic equations used to describe self-gravitating fluids with a mag- 
netic fleld. Throughout this paper, we assume perfect flux freezing. Also, we discuss an 
equation of state of particular interest — a non-isothermal barotropic equation of state that 
heuristically includes the effects of the turbulence that is observed in molecular clouds. 

In dimensionless units (see Appendix ^), the equations of ideal fluid dynamics with 
self-gravity and a frozen magnetic fleld can be written in the form 

|^ + V-(pu) = 0, (1) 

p— ^ + p(u- V)u + Vp + pV^ - (V X B) X B = 0, (2) 

ot 

^ + V X (B X u) = 0, (3) 

= p, (4) 

where p is the mass density, u is the fluid velocity, p is the pressure, ip is the gravitational 
potential, and B is the magnetic fleld strength. We will assume the pressure to have a 
barotropic form, i.e., 

p = P{p). (5) 
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Most previous theoretical studies of molecular clouds have assumed an isothermal equa- 
tion of state, p = c^p, where Cg is the thermal sound speed. We generalize this equation of 
state by adding a term which attempts to model the "turbulence" observed in molecular 
clouds. The equation of state takes the form 

p = clp + po\og{p/p), (6) 

where po is a constant, which may be determined empirically, and p is an arbitrary reference 
density ([Lizano fc Shu 1989|) . In this paper, we will always take our dimensionless density 



variable p to be 1 at the center of equilibrium configurations, so that p = pc, the central 
density (see Appendix 0). The logarithmic nature of the turbulence term arises from the 
empirical linewidth-density relation Aw oc p~^/^ ( [Larson 1981| ; [Myers 1983| ; [Dame et al. 



[1986[ ; [Myers fc Fuller 1992| ). After transforming to dimensionless variables, this equation of 
state can be written as 

p = P(p) = p + K\ogp, (7) 
where the turbulence parameter n is defined by 

f^ = Po/clp. (8) 

Typically, we are interested in clouds with number densities of about 1000 cm~'^ and thermal 
sound speed Cg ~ 0.20 km sec~^. Observational considerations suggest that po ranges from 
10 to 70 picodyne/cm^ ([Myers fc Goodman 1988 ; Solomon et al. 1987 ). Thus, molecular 



clouds are expected to have values of the turbulence parameter in the range 6 < k < 50. 
In this work, we want to include both the isothermal limit and the pure "logatropic" limit 
where the turbulence dominates, so we use the expanded range < k < oo. 

In addition to the empirical motivation for using an equation of state that is softer 
than isothermal, there exists theoretical support as well. It has been shown that when a 
molecular cloud begins to collapse (on large spatial scales), a wide spectrum of small scale 
wave motions can be excited ( [Arons fc Max 1975[ ; see also [Elmegreen 1990| ). Additional 
energy input and wave excitation can also be produced by outflows from forming stars 
( [Norman fc Silk 1980[) . In the presence of magnetic fields, these wave motions generally take 



the form of magnetoacoustic and Alfven waves. In any case, these small scale wave motions 
have velocity perturbations which vary with the gas density in rough agreement with the 
observations described above. These wave motions can be modeled with an effective equation 
of state which is softer than isothermal (see [batuzzo fc Adams 1993[; [McKee fc Zweibel 1995[). 



- 5 - 



3. Static Equilibrium 

We take the unperturbed state of the system to be one of hydrostatic equihbrium 
(u = 0). Since the unperturbed magnetic field is uniform, it provides no support. There- 
fore, the unperturbed state is unaffected by the presence of the magnetic field and, hence, 
is identical to that discussed in [Paper I| . The equilibrium density distribution is determined 
by the equation 

-V^p - \vp ■ Vp + p = 0. (9) 
P P^ 



3.1. Filaments 

For the filament, we adopt cylindrical coordinates (r, 0, z) and assume azimuthal sym- 
metry. The equilibrium equation becomes 



d'^p 1 dp 
dr'^ r dr 



P"{P) 1 
P'(p) p. 




0, (10) 



where the primes denote derivatives with respect to density, e.g. -P'(p) = dp/ dp. In the 
isothermal case, p = P{p) = p, the equilibrium solution is well known ( pstriker 1964| ) and 
has the simple form 

p(r) = (l +rV8)"'. (11) 

For equations of state which include turbulence (k ^ 0), the solutions have been found 
numerically (see Paper I| ). 



Paper 1 demonstrated an important difference in the asymptotic behavior of the equi- 
libria of the isothermal and turbulent cases. At large radii, the equilibrium density profile 
of the isothermal filament behaves as r~'^ while that of the turbulent filament falls off only 
as r~^. This result implies that the mass per unit length of the isothermal filament remains 
finite, whereas that of the turbulent filament diverges. The mass per unit length of the 
isothermal filament is given in dimensionless units by p = Svr. In addition to providing an 
observable test of the significance of turbulence in clouds, these differences have important 
implications for the mass scales of fragmentation. 



3.2. Sheets 

For the thick sheet (sometimes referred to as a slab or stratified layer), we use cartesian 
coordinates {x, y, z) and assume translational symmetry in the y- and z-directions. The 
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equilibrium equation now becomes 



+ 



P"{p) 
P'ip) 




P'ip) 



0. 



(12) 



(13) 

and for turbulent equations of state the solutions have been found numerically ( Paper l|) . 
As in the case of filaments, the asymptotic behavior of the isothermal and turbulent sheets 
differ significantly. [Paper 1] showed that the surface density of the isothermal sheet is 2^/^ 
and that of the turbulent sheet is infinite. 



Again, the isothermal solution is well known ( |Spitzer 1942| ; |Shu 199^ ) 

p = sech^ (^x 



4. Perturbations 

We now consider the addition of small perturbations to the equilibrium states described 
in We denote the equilibrium quantities with a '0' subscript and the perturbations with 
a '1' subscript. After linearization, the fluid equations (0^ can be written 

dpi 



dt 



+ V ■ (poui) = 0, (14) 



Po-^ + + PoV^i + piV^o - = 0, (15) 
at 

^ + V X (Bo X ui) = 0, (16) 

V'V'i = Pi, (17) 

where £i is the first-order mean Lorentz force per unit volume acting on the fluid and is 
given by 

Ci = (Bo ■ V)Bi + (Bi ■ V)Bo - V(Bo ■ Bi). (18) 
Adopting a barotropic equation of state (|^), we may write the pressure perturbation as 

Pi = P'iPo)pi. (19) 
We assume that the equilibrium magnetic field is axial, i.e.. 

Bo = 5o(^)z, (20) 
and seek solutions for the perturbations of the form 

pi(x, t) = /(tu) exp(i/c2; — zcut), (21) 
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Ui(x, t) — v(ro) exp{ikz — icut), 
Bi(x, t) = h{zu) exp{ikz — iut), 
■i/'i(x, t) = (f){zu) exp{ikz — icut), 



(22) 
(23) 

(24) 



where cc7 is a generalized coordinate, which we take as our cartesian x in the case of the 
sheet and our cyhndrical r in the case of the filament. We find purely algebraic expressions 
for 6^, bz, and v/. 



k ^ 
iBo ( 1 dpo 



9 u ; J 9 



1 dBo k^BodBo\ 
\- — ; — I v^, 



^Po dva Bq duj uj^Po d'^ 
iBq dBo \ 



(25) 
(26) 

(27) 



u \po ujpo dw j 

and eliminate them from the remaining differential equations. Following a convenient change 
of one variable, 

w = iuJVj^j, (28) 
we obtain the following ordinary differential equations: 



pVw + {uj' - k'P')f - k'p(t) + 



'dp e dB\ 

dw uj^ dw I 



(29) 



dw 

V-^ - k\ 

dw 



di 



p - —B' -7- + ^20 + = 0, 
I o;^ / dw 



/ = 0, 



(30) 



(31) 



where we have omitted the '0' subscripts on equilibrium quantities. We have introduced the 
general differential operator 



V 



the coefficient functions 

k^P'B^dB ^B^ 




A, 



uj'^p 



oj^p^ dw 

2B^\ dp dip 



cartesian 



cylindrical 



^ + ^]^+ (p" _ ^' 
p B j dw \ p ) dw 



+ P" 



+ 



3BdB 



dw dw p dw ' 



(32) 



(33) 
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A, 



A, 



k^B^dB e 



cu^p dzu 



k^B^ 
+ — 




3 dB' 
B dzu 



(34) 




k^B^dB fd^B IdB 4 rfB 
tu^p dw I dcc^ / (itu B dw 




k' + - 



1 rfV 



p (itz?^ 



3 d p dB 



C 



pB dw dw \ B dw 
and the cylindrical coefficient function 



\ p dw J B dw"^ 
1 dB 



1 dB 
B dw 




C = V 



d 

dw 



For boundary conditions, we take 

1 ^ = 
dw 

f = 0,w 



f 



cartesian 



- cylindrical 

r 



0, w = at Ci7 = 
= at = 00. 



(35) 



(36) 



(37) 



Real molecular clouds will, of course, be pressure confined at some finite distance from the 
center. However, as long as this distance is larger than the scale height (in the w direction), 
the influence of the pressure confinement will be small (see [Nagasawa 1987|). 



5. Perturbation Solutions and Dispersion Relations 

In this section, we discuss the solutions and dispersion relations for the perturbation 
problem presented in the preceding section. Equations (p9|-|3TD constitute a disguised eigen- 
value problem; the equations could be cast into the usual eigenvalue equation form. For 
numerical purposes, we choose, instead, a form which isolates terms involving the deriva- 
tives of the perturbation. Also, in this form, the possibility of a singularity becomes evident. 
The coefficient of d f / dw vanishes when 

p - t^—.B^ = 0, (38) 

where we have adopted the equation of state (^. In the isothermal case, a singularity occurs 
where 

P-B'g-l). (39) 



- 9 - 



This singularity indicates that the addition of a magnetic field can preclude the existence 
of propagating (c<j^ > 0) wave solutions that are localized in the filament. This absence of 
propagating waves occurs because Alfven waves on the uniform magnetic field background 
can easily carry the energy of a propagating wave away from the filament, causing the wave 
to decay in time. We avoid this complication by solving the problem only in the unstable 
regime, where iJ^ < 0, and hence no singularity arises. 



For the marginally stable perturbation, equations p9|-|3lD can be simplified by applying 
the conditions c^;^ = and w = 0. When the equilibrium magnetic field strength B is 
uniform, we find 

^'7o + P0o = O, (40) 

- > - («) 

We solve this well known equation to determine the critical wavenumber and associated 
marginally stable perturbation functions, /o and (po. 

Once we have the marginally stable solution, we proceed to smaller values for the 
wavenumber k. We use a relaxation algorithm with a finite difference approximation ( [Press 



et al. 199"^ ) to solve equations (p9|-pl|) for successively smaller values of k , where the solution 



at the previous step is used as an initial approximation. 



5.1. Filaments 

First, we consider the case of molecular cloud filaments. We adopt cylindrical coordi- 
nates (r, 0, z) and assume azimuthal symmetry. We use the one-parameter turbulent equation 
of state (^. Notice that the limit of vanishing turbulence parameter k — corresponds to a 
purely isothermal cloud. The opposite limit, k ^ oo, corresponds to a highly turbulent cloud 
with a purely logatropic equation of state p = log p (see Appendix ^ for details of the differ- 
ences in the scaling of variables for the purely logatropic case). Observational considerations 
indicate that the range 6 < k < 50 is the most physically relevant. We also assume that the 
unperturbed magnetic field is uniform. Using typical values for the observed characteristics 
of molecular clouds, B ~ lOyuG, p ~ 4 ■ 10^^° g cm~^, and Cg ~ 0.20 km sec~^, we find that 
the physically interesting magnitude of the dimensionless field strength variable is i? ~ 0.7. 

All the dispersion relations possess the same general form and exhibit a minimum value 
of u"^ at a certain wavenumber, which we will call fcfast- Since u"^ < throughout the region 
that we explore, all solutions that we find represent unstable modes. The magnitude of 
the frequency \uj\ is the rate of growth for a particular mode, and hence the minimum of 
the dispersion relation corresponds to the fastest growing mode. The location of this point 
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(^fastj ^^fast) dispersion plane varies as the parameters k and B change. In Figure |l|a, 

we show isothermal (k = 0) dispersion relations for various values of the equilibrium field 
strength B. These results agree well with corresponding results of Nagasawa (1987), which 
were obtained using a different numerical technique. Notice that the magnetic field tends to 
decrease the instability of all unstable modes but has no influence on the critical wavenumber 
and has only a very small effect on the fast wavenumber fcfast- This effect is contrary to that 
exhibited by models which adopt a non-uniform unperturbed magnetic field. Nakamura et 
al. (1993) show that an isothermal filament with a constant ratio of magnetic pressure to 
thermal pressure is destabilized by increasing the strength of the magnetic field. 

Figure ^ shows dispersion relations for various values of the turbulence parameter k 
with no magnetic field. Notice that the wavenumber has been scaled by the factor (1 + k)^/^, 
which is just the effective sound speed dP/dp at the center of the cloud. This scaling is 
done to account for the fact that the variable transformation performed to eliminate units 
uses the thermal sound speed. Displayed this way, the dispersion relations can be seen to 
converge to the limiting case where p = log p as the turbulence parameter k, increases. Again, 
these results agree well with corresponding results in [Paper 1| . Contrary to the effect of the 
magnetic field, the presence of turbulence significantly increases the growth rate of the fast 
mode as well as decreases its scaled wavenumber. In addition, turbulence is seen to decrease 
the scaled critical wavenumber. [Paper 1| closely examines these effects for turbulent clouds 



without magnetic fields. 

Figures 0b and similarly show more dispersion relations where one parameter is 
varied while the other is held fixed. In these figures, however, the fixed parameter (the field 
strength B or the turbulence parameter k) has a more physically relevant value. Notice in 
Figure 0b that the magnetic field is not as effective at decreasing the instability as in the 
isothermal case (Figure 0a). This result can be understood by again considering the effective 
sound speed. The dimensionless magnetic field strength variable scales inversely with the 
thermal sound speed. A more analytically convenient variable would use the effective sound 
speed instead of the thermal sound speed. We choose to use the thermal sound speed in 
our transformation because of its physical significance; the thermal sound speed is well 
determined in observations of molecular clouds. The decreased effectiveness of the magnetic 
field is shown clearly in Figure ^, which plots the growth rate of the fastest growing mode 
versus the equilibrium field strength for a variety of values of the turbulence parameter. For 
the isothermal filament, the effect of the magnetic field saturates around a field strength of 
i? ~ 2, whereas a field strength of -B 10 is required to achieve the maximum effect for the 
turbulent filament with k = 10. 



The wavenumber of the fastest growing mode defines a length scale Afast = '^T^/hast 



- 11 - 



for the fragmentation of the filament. This physical length scale is tabulated in Table |l| for 
several physically interesting values of the magnetic field strength and turbulence parameter. 
Notice that the influence of turbulence dominates the effect of the magnetic field. For 
isothermal filaments, a corresponding mass scale of fragmentation is set by 

Mfrag = STrAfast. (42) 

Paper 1| showed that the physical mass scale for the isothermal filament is given by 

Pc 



Mfrag ~ 14.5 Mq 



0.20 km/s 



4 X 10-20g/cm^ 



-1/2 



(43) 



Here, in the strong magnetic field limit, we find this mass scale to be slightly decreased; 

3 r H -1/2 



frag 



13.5 Mr7 



© 



Pc 



(44) 

0.20 km/s ^ on , 3 ■ \ I 

Since the mass per unit length of turbulent filaments is formally infinite, so is their fragmen- 
tation mass scale. 

Measuring the actual physical spacing of clumps along a molecular cloud is rather dif- 
ficult, since the uncertainty in the distance and orientation of the cloud is often large. A 
more easily measured quantity is the aspect ratio 

T= , (45) 

-Rhwhm 

where -Rhwhm is the half-width at half- maximum of the filament. The dependency of this 
ratio on the parameters k and B is illustrated in Figure ^. Notice that the effect of the 
magnetic field on the aspect ratio JF is rather small for isothermal clouds compared to the 
effect on turbulent clouds. We also note that there is a threshold field strength, below which 
there is little change in JF and above which the influence of the magnetic field increases 
dramatically. 

Figures |^ and ^ show cross-sections of the perturbed density and the flow of the fluid 
and the magnetic field lines for various filaments. Isothermal filaments with and without 
magnetic field are shown in Figure ^; similarly, turbulent filaments are shown in Figure ^ 
These figures illustrate the physical reason for the limited effectiveness of the magnetic field. 
For B = 1, the effect of the magnetic field in the isothermal filament is near saturation. This 
saturation is manifested by the nearly axial flow of the fluid. In contrast, when the magnetic 
field is absent the flow has a significant radial component. Notice that this contrast is not 
present in the turbulent filament, where the effect of the magnetic field at S = 1 is still far 
from saturation. These results support the conclusion that the saturation occurs because 
the magnetic field cannot prevent contraction of the fluid along the field lines ( [Nagasawa 
19871) . 
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5.2. Sheets 

Now, we consider a molecular cloud sheet or stratified layer. We adopt cartesian co- 
ordinates {x, y, z) and assume translational symmetry in the y-direction. Here, vj = x, 
the differential operator V is just a derivative, and the cylindrical coefficient function C is 
identically zero. The problem is solved using the same numerical technique as used for the 
filament. 

The dispersion relations for the sheet (Figures |^ and ^ possess the same general form 
and exhibit similar dependencies on the magnetic field and turbulence parameters as those 
for the filament. In general, the growth rates are about 30-40% larger for the sheet than 
for the filament. The infiuence of the magnetic field and turbulence on the fastest growth 
rate is displayed in Figure ||. Notice that these effects are qualitatively the same as for the 
filament discussed in the previous section. The magnetic field tends to stabilize the cloud. 
As before with filaments, this effect is contrary to that exhibited in isothermal sheets with a 
constant ratio of magnetic pressure to thermal pressure ( [Nakano 1988| ; [Nakamura, Hanawa 
fc Nakano 19911) . 



We tabulate the physical fragmentation length scale Afast in Table 0. As before, for the 
isothermal filament, this length scale yields a mass scale of fragmentation for the isothermal 
sheet. The magnetic field, however, does not significantly affect the the length scale of the 



isothermal sheet and, hence, has little effect upon the corresponding mass scale (see Paper I 
for a discussion of the mass scales) . 



6. Discussion &: Summary 

In this paper, we have studied gravitational instabilities of astrophysical ffuids in two 
spatial dimensions. Although many of the results are general, molecular clouds provide our 
primary motivation. We have accounted for the presence of turbulence in molecular clouds 
by adding a logarithmic term in the equation of state and included the inffuence of a uniform 
axial magnetic field perfectly frozen in the fiuid. We have studied wave-like perturbations on 
the hydrostatic equilibrium states and have numerically determined the dispersion relations 
and the structure of the perturbations in the unstable regime. 
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6.1. Summary of Results 

We have found a number of results concerning gravitational instabilities in molecular 
clouds, which add to our general understanding of fluid dynamics in self-gravitating systems 
with uniform magnetic flelds. Our main results can be summarized as follows: 

1. The dispersion relations for these two dimensional perturbations have the same general 
form for all values of magnetic fleld strength and for all equations of state considered 
here. In particular, the growth rate for gravitational instability obtains a maximum 
value for one particular wavelength of the perturbation (see Figures and 0). 

2. Turbulence increases (~ 22%) the growth rate of the fastest growing mode when the 
magnetic fleld is not strong (see Figures || and |). 

3. Turbulence increases the length scale of the fastest growing mode when the magnetic 
fleld is not strong. Furthermore, turbulence increases the aspect ratio T of the structure 
by about 29% (see Figures || and H). 

4. A magnetic fleld tends to decrease the growth rate of all unstable perturbations (see 
Figures ^ and The growth rate of the fastest growing mode is suppressed by 
about 12% in the absence of turbulence. This suppression is much stronger (~ 31%) 
in turbulent clouds. Since the fastest growth rate is signiflcantly larger, however, 
for turbulent clouds without magnetic flelds, the growth rate is nearly the same for 
isothermal and turbulent clouds in the strong field limit. 

5. A magnetic fleld has a relatively small influence (~ 6%) on the length scale of fragmen- 
tation for isothermal clouds. For turbulent clouds, there is httle effect on the length 
scale unless the magnetic fleld exceeds a threshold, which is of the order of the turbu- 
lence parameter. Above this threshold, the magnetic fleld reduces the fragmentation 
length scale of turbulent clouds up to about 15% (see Figure]^). 

6. In certain regimes, turbulence and magnetic fleld effects counteract one another, in 
the sense that turbulence enhances the growth while the magnetic fleld suppresses the 
growth. For example, a molecular cloud fllament with k = 10 and B = 2 has essentially 
the same fragmentation rate |a;fast| as a fllament with k = = B (see Figure 

7. The mass scale for fragmentation of isothermal fllaments depends only weakly on the 
magnetic fleld strength and typically has a value of ~ 14 Mq, which is much larger 
than a typical stellar mass. The mass scale for a turbulent fllament is formally inflnite 
(see also [Paper 1| ). 
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8. The dispersion relations for molecular cloud sheets are qualitatively similar to those 
of filamentary clouds. Moreover, the fragmentation rate \ujiast\ and length scale Afast 
depend on the turbulence parameter k and the magnetic field strength B in the same 
way as in the case of the filament (compare Figures ^ and ^). 



6.2. Implications for Observed Molecular Clouds 

The results of this paper can be compared with observed structures in molecular clouds. 
In particular, these theoretical calculations make several definite predictions that can be used 
to further our understanding of the physics of molecular cloud systems. These predictions 
are outlined below. 

For given values of the turbulence parameter k and the magnetic field strength B, the 
dispersion relations predict that gravitational fragmentation occurs with a well defined length 
scale Afast and has a well defined growth rate |ci;fast|- The length scale Afast can be deduced 
from observed molecular clouds maps and compared with these theoretical results. If the 
turbulence parameter k can be independently determined from observations of molecular 
line-widths, and if the magnetic field strength B can also be measured, then one can obtain 
a direct test of the theory. 

Unfortunately, however, observational uncertainties often prevent a clean determination 
of the length scale Afast- These uncertainties arise from the (often unknown) projection angle 
and the (imprecisely known) distance to the cloud. This latter effect can be avoided by 
considering the aspect ratio JF = Afast /-Rhwhm of the perturbation. This ratio is shown as a 
function of magnetic field strength B and the turbulence parameter k, in Figure ^. 

Observations indicate that the strength of the magnetic field in actual molecular clouds 
increases with increasing mass density ( [Troland fc Heiles 1986| ; [Heiles et al. 1993| ), rather 



than remaining constant throughout the cloud as we have assumed in our model. Since the 
equilibrium configuration of such clouds are supported in part by the magnetic field, this 
difference can lead to qualitative changes in the dependence of the fast mode on the strength 
of the field. For example, isothermal models that adopt a constant ratio of magnetic to 
thermal pressures, i.e., B oc p^^^, show that the magnetic field destabilizes the cloud and 
decreases the length scale of fragmentation ([Nakamura, Hanawa, fc Nakano 1993| ; [Nakamura 



[Hanawa, fc Nakano 19911) . Several reasonably successful efforts to apply such models to 



particular clouds have been made ( [Hanawa et al. 1993 ; Nfakamura, Hanawa, fc Nakano 19"93 



[Matsumoto, Nakamura, &: Hanawa 1994|) . 



Another important observational diagnostic is the structure of the velocity field. 



Paper 1 
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showed that unstable (i.e., growing) perturbations have a completely different velocity field 
signature than propagating waves. Furthermore, as we show in this paper, the velocity 
field of the perturbations depends on both the magnetic field strength and the turbulence 
parameter (see Figures ^ and In particular, as the magnetic field strength increases, the 
radial velocities become smaller and the flow is almost entirely in the z direction (along the 
filament - see Figure |^). Observations of the geometry of the velocity field in molecular 
cloud filaments can thus be a powerful diagnostic of the internal dynamics. 



6.3. Directions for Future Work 

In this paper, we have begun to study the evolution of structure in molecular clouds. 
Numerous directions for future work, however, remain to be explored. These directions 
include both theoretical and observational studies. 

Small scale magnetohydrodynamic motions, such as Alfven waves, are often considered 
to account for the turbulence observed in molecular clouds. Here, for the first time, we have 
attempted to study the effects of this turbulence in conjunction with the effects of large 
scale magnetic fields, which are widely believed to provide a significant fraction of pressure 
support for molecular clouds. Our model, however, assumed a uniform unperturbed magnetic 
field, which provides no support to the equilibrium state of the cloud but does impede 
the fragmentation of the cloud by gravitational instability. A more relevant unperturbed 
configuration should include a magnetic field which does contribute to the support of the 
equilibrium state. Observations of molecular clouds suggest that the field strength often 
varies with the density roughly according to i? ~ p^^^ (e.g., see [I'roland fc Heiles 1986| and 
Heiles et al. 1993| ). Furthermore, this relation is theoretically attractive since it implies a 



constant ratio of magnetic to thermal energy. Stodolkiewicz (1963) determined the critical 
length scale (27r/A;o) for such an isothermal system. Though several studies on isothermal 
clouds with such an unperturbed magnetic field have been performed ([Nakamura, Hanawa, fc 



Nakano 1993| ; [Nakamura, Hanawa, fc Nakano 1991|) , the effects of turbulence in the presence 



of such a field have yet to be investigated. 

In future work, it will be interesting to investigate the behavior of perturbations using a 
wider variety of equations of state. We have focused here on equations of state that are softer 
than isothermal; stiffer equations of state should also be examined. For these equations of 
state, the equilibrium configurations of filaments have a finite radius, outside of which the 
density is zero. Other functional forms for the equation of state of a turbulent gas should 
be considered as well. The logarithmic term used in this work was deduced empirically. 
Theoretical studies of turbulence in molecular clouds may lead to different formulations for 
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the turbulent contribution to the pressure. 

Finally, we note that the clumps observed in molecular cloud filaments are generally not 
small perturbations on an equilibrium filament. In general, molecular cloud substructure 
lies in the fully nonlinear regime. In previous work, we began to study nonlinear waves in 
molecular clouds (e.g., [Adams, Fatuzzo, fc Watkins 1994| ; see also [Infeld fc Rowlands 19901 ). 



This work was restricted to one spatial dimension, although two-dimensional effects were 
heuristically taken into account. Future work should study two-dimensional wave motions 
and instabilities in the fully nonlinear regime. 

We would like to thank Greg Laughlin, Marco Fatuzzo, Lee Hartmann, Phil Myers, 
Paul Ho, and Jennifer Wiseman for constructive comments and discussion. This work was 
supported by a NSF Young Investigator Award, NASA Grant No. NAG 5-2869, and by funds 
from the Physics Department at the University of Michigan. 



A. Transformation to Dimensionless Variables 

The equations of motion of the physical fields of a fluid with a perfectly frozen magnetic 
field can be written as 

|^ + V-(pu) = 0, (Al) 
p-^ + p(u ■ V)u + Vp + pVtp - —(V X B) X B = 0, (A2) 

^ + V X (B X u) = 0, (A3) 

= 47rG'p, (A4) 

where p, u, p, and ip are the mass density, velocity, pressure, and gravitational potential of 
the fluid, respectively, and B is the magnetic field. To simplify the problem, we transform 
to dimensionless variables according to 

t —>■ tt p ^ pp u ^ iiu 

X ^ t-ux p ^ pv?p B ^ (47rp)^/^{iB (A5) 
ip — > V?1p 

where i = (47rGp)^^/^. Throughout this paper, we let p = pc, the central density of the 
equilibrium state. For cases where the equation of state takes the form of equation (0), 
we let u = Cg, the thermal sound speed. When we consider the purely logatropic equation 
of state p = p\og{p/p), we set u = {p/pY^"^- Under the transformation (|A5|) , the fluid 
equations (|A1|-|A^) are cast into the form (|1|-^. 
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Table 1. Physical Fragmentation Length Scale for Filament 



Magnetic Field 


Turbulence parameter 


K 




B 


1 


5 


10 


20 





0.777 1.51 


2.89 


3.99 


5.57 


0.5 


0.752 1.52 


2.93 


4.02 


5.59 


1.0 


0.738 1.42 


2.96 


4.09 


5.65 


2.0 


0.731 1.31 


2.74 


4.01 


5.74 


10 


0.729 1.26 


2.39 


3.32 


4.69 


Note. — The values in this table give the fragmentation length scale 


in parsecs. 


The values scale with the thermal sound speed and 


mass density according to A ~ 


(cs/0.20 km sec 


^) (pc/4.0 X 10-20 g cm-3) 










Table 2. Physical Fragmentation Length Scale for Sheet 




Magnetic Field 


Turbulence parameter k 






B 


1 


5 


10 


20 





0.677 1.10 


2.02 


2.76 


3.84 


0.5 


0.677 1.12 


2.03 


2.78 


3.85 


1.0 


0.677 1.13 


2.07 


2.81 


3.87 


2.0 


0.675 1.12 


2.10 


2.87 


3.95 


10 


0.673 1.08 


1.99 


2.75 


3.87 



Note. — The values in this table give the fragmentation length scale in parsecs. 
The values scale with the thermal sound speed and mass density according to A ~ 
(cs/0.20 km sec-i) (pc/4.0 x 10-20 g cm-^)"^/^ 
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Figure Captions 



Fig. 1. — Dispersion relations for molecular cloud filaments. The abscissa is the wave number 
multiplied by the effective sound speed, and the ordinate is the square of the frequency. For 
any particular wave number in the region plotted, the growth rate (|a;|) decreases with 
increasing magnetic field strength, (a) Isothermal filament [k = 0) with magnetic field 
strengths of 0, 0.05, 0.10, 0.20, 0.50, 1.0, and 2.0. (b) Turbulent filament {k = 10) with 
magnetic field strengths of 0, 0.20, 0.50, 1.0, 2.0, 5.0, and 10. 

Fig. 2. — Dispersion relations for molecular cloud filaments. Solid curves show isothermal 
dispersion relations. The abscissa is the wave number multiplied by the effective sound 
speed at the center of the filament. The ordinate is the square of the frequency. Dispersion 
relations are shown for various values of the turbulence parameter k. Notice that the location 
of the minimum moves to larger growth rates and smaller wavenumbers as the turbulence 
parameter k, increases, (a) Filament with no magnetic field. «; = 0, 1, 2, 5, 10, oo. (b) 
Magnetic filament {B = 1). k = 0, 1, 2, 5, 10, 20. 

Fig. 3. — Growth rate of the fastest growing mode. Points represent minima of dispersion 
relations for various magnetic field and turbulence strengths. Points for fixed values of the 
turbulence parameter {k = 0, 1, 2, 5, 10, 20) are interpolated with cubic splines. 

Fig. 4. — Relative length scale of fragmentation. Points represent minima of dispersion 
relations for various magnetic field and turbulence strengths. The ordinate is the aspect 
ratio, the wavelength of the fastest growing mode divided by the half-width at half-maximum 
of the unperturbed filament. Points for fixed values of the turbulence parameter [k — 0, 1, 
2, 5, 10, 20) are interpolated with cubic splines. 

Fig. 5. — Cross-section of the fastest growing mode in isothermal filaments without and with 
magnetic field {B — 0, 1). The total density (equilibrium plus perturbation) is represented 
with solid contours (p = 0.2, 0.4, 0.6, 0.8, 1.0, 1.2), and the fluid velocity is represented by 
the array of arrows on the right half of the filament. When present, magnetic field lines are 
shown as dashed lines on the left side. 

Fig. 6. — Cross-section of the fastest growing mode in turbulent filaments {k, — 10) without 
and with magnetic field {B = 0, 1). The total density (equilibrium plus perturbation) is 
represented with contours (p = 0.2, 0.4, 0.6, 0.8, 1.0, 1.2), and the fluid velocity is represented 
by the array of arrows on the right half of the filament. When present, magnetic field lines 
are shown as dashed lines on the left side. 
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Fig. 7. — Dispersion relations for turbulent molecular cloud sheets with magnetic field 
strengths of 0, 0.5, 1.0, 2.0, 5.0, and 10. The turbulence parameter k = 10 is fixed. For 
any particular wave number, the growth rate decreases with increasing magnetic field 
strength. 

Fig. 8. — Dispersion relations for molecular cloud sheets with a fixed magnetic field strength 
oi B — 1.0 and various degrees of turbulence {k, — 0, 1, 2, 5, 10, 20). The abscissa is the wave 
number multiplied by the effective sound speed at the center of the filament. The fastest 
growth rate (|a;fast|) increases with increasing values of the turbulence parameter k. 

Fig. 9. — Growth rate of the fastest growing mode in the sheet. Points represent minima 
of dispersion relations for various magnetic field and turbulence strengths. Points for fixed 
values of the turbulence parameter {k — 0, 1, 2, 5, 10, 20) are interpolated with cubic splines. 
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